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1. Introduction 

The study of the traversal of potential barriers by wave packets in the 1-dimensional case 
is relevant for the electron transport through barrier junctions, the physics of wave guides 
and light transmission through Photonic Band Gaps, etc. On more theoretical grounds, 
the resonances, transient excitations and the time of arrival of wave packets, with ensuing 
paradoxes of non-locality and super-luminal tunnelling, are further interesting aspects. 

Relevant resonances are mostly expected in systems of two or more barriers 
separated by a gap. Actually they generally occur in any simple barrier, and the plain 
square barrier is most tractable, fully representative for many theoretical purposes and 
devoid of bound states or other unessential features. 

On the other hand, practical calculations face important limitations and the effort 
of pushing the analytical calculation as far as possible is always rewarding. Significant 
progress has been accomplished [1] for a Gaussian packet impinging on a square barrier, 
where the Gaussian structure was exploited, and [2] for the shutter problem, where the 
contribution of the resonant (Gamow) states in a double square barrier was worked out. 

In this work we consider the time evolution, in the presence of a square barrier 
of height V in the interval [0, L], of an initial plane wave with support on the whole 
1-D space, aiming to work out the contribution of the resonant states. Our calculation 
extends the results of the shutter problem, where the initial plane wave occupied the 
space axis at one side of the barrier and was specially suited to study details of the 
propagation of the wave front. 

The reason for doing this is twofold. First, the simpler problem of the shutter initial 
condition ip s (x) = 6(—x)e lkx has been assumed in the literature to be meaningful also 
for the analysis of the transmission in the tunnelling regime when k -C ky = \/2mV . 
Unfortunately such a cut-off initial wave function is highly non-monochromatic as its 
Fourier transform shows: 

Mp) = ^=T dxe-**rl> a (x) = + t/| S(k - p) 

V 271 J —oo 1 y 271 K — p V 2 

The Principal Part distribution around k is exceedingly wide so that the amplitude 
of the momentum p = ky skimming the barrier is oc (k — ky)~ x , a sizeable value (in 
modulus, actually it is > ky 1 ) even for very small incident k , with a slow decrease for 
p > ky . Secondly, by considering infinite plane waves the way is paved for further uses, 
since superpositions of them build up any desired initial wave packet. In particular, the 
compact result obtained for the Gaussian wave packet P could be re-derived pinpointing 
the contribution of each single Fourier component. 

In this work we follow the approach used by Peierls and Garcia-Calderon [3], 
namely Laplace-transforming the time evolution equation into a second order linear 
differential equation (SOLDE) in one variable, expressing its solutions in terms of the 
Green function (GF) with resonant boundary conditions (RBCs) and then undoing the 
Laplace transform. 
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The first step is done in Section 2, where a brief discussion is made of the use of 
the Green's method for the solution of SOLDEs when the boundary conditions (BCs) 
required for the GF and for the solution are different. The GF with RBCs features a 
simple structure of isolated (resonance) poles, in terms of which the (p-dependent) 
Laplace-transformed wave function is worked out in Section 3. Then the explicit 
calculation of the inverse Laplace transform leading to the sought after t-dependent 
solution is carried out in Section 4. 

As a check of this cumbersome analytical calculation, the t — (Section 5) and the 
t = oo (Section 6) limits of the time-dependent solution i/)(x,t) are calculated: at t — 
one must recover the initial plane wave and at t = oo a suitable stationary solution 
must be reached. Then the conclusions are drawn in Section 7. Some notations and a 
number of calculations and technical details are deferred to the Appendices. 



2. The evolution equations 

We consider the 1-D time- dependent Schrodinger equation 

(ih^-H)rl>(x,t) = (1) 

where the Hamiltonian H = p 2 / 2m + V(x) corresponds to the square barrier potential 
V(x) = 9(x)9(L — x)V and the solution satisfies the initial condition ip(x, 0) = ip (x) = 
e ikx , k > . 

The Laplace transform ip(s) = J °° dt e~ st ip(t) on the time variable, brings the 
parabolic partial derivative differential equation (1) to the simpler SOLDE 

<9 2 2m 

L x $(x, p) = i—+p 2 - -jfV(x)} p) = ia e ikx (2) 

where a = 2m /h and p 2 = ia s (Appendix A). Besides the usual scattering in and out 
solutions, the time-independent Schrodinger homogeneous equation L x ijj = has the 
resonant solutions satisfying the homogeneous outgoing RBCs 

d x i)\ x =o= -ip^(O) , d x i[)\ x=L = ipi>(L) , (3) 

which exist only for a denumerable set of isolated values p n of p lying in the lower half 
complex plane (Appendix A). 

A solution to the inhomogeneous equation (2) for p 2 e R+ can be explicitly written 
down for each region of V(x) : 

Mx,p) =Be~ ipx + ^—e ikx (x<0) 

P rC 

$n(x,p) = M e ip ' x + N e~ ip ' x + - ^ , 9 e ikx (0 < x < L) (4) 

$ IH (x, p) = A e ipx + -^^^ (x > L) 

where the amplitudes A, B, M and iV are functions of p completely determined by the 
matching conditions at x = and x = L and have a common denominator of the form 
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D{p) = Q 2 e ip ' L — (B 2 e~ ip ' L (notation in Appendix A). Out of all the terms of the general 
homogeneous solution, the choice of e ipx (for x > L) or e~ ipx (for x < 0) is dictated by 
the behaviour of the solution for x — > ±00 , since a small positive imaginary part in p 
is supposed when performing the inverse Laplace transform leading to ip(x,t) , t > 
(Appendix D). 

However, this expression of the solution is not suited to perform the inverse Laplace 
transform back to tp(x,t) because of the nontrivial analytical form of the amplitudes 
above. Instead the GF approach lets us to express 4>(x,p) as a sum of isolated pole 
terms plus other simple terms easier to deal with. To this end, (2) together with the 
Green equation for G(x,x',p) can be written in the form 

L x $(x,p) -iae ikx =0 (5) 

L x G(x,x',p)-S(x-x') =0 , (6) 

where G(x,x',p) is required to obey RBCs, namely 

d x G(x,x',p)\ x=0 = -ipG(0,x',p) , d x G(x, x',p)\ x=L = \pG(L, x',p) (7) 

It is expected that G(x, x',p) = L^ 1 will have poles at p = p n , where the homogeneous 
equation L x ip = with RBCs has non-trivial solutions. 

The equations (5) and (6) may conveniently be given the short-hand notation 
^ = and T = respectively. Then the Green method uses the integral equation 
J L dx [ipT — G^} = to obtain vp(x,p) in terms of G(x, x',p). Only if ip and G obey the 
same homogeneous BCs, the surface terms in this integral cancel out and one obtains 
the familiar result ipnix, p) = J L dx' ia G(x', x, p) e lkx ' . 

However the BCs for ip(x,p) (involving ip and d x ip both at x = and at x — L) 
stemming from (4) are different from (3) and non- homogeneous. In that case one obtains 
a modified expression for the region //, namely 

\kx' 



ipn(x,p) = / dx' ia G(x' , x,p) e 1 ' 
Jo 



a _ ifeL a 



k G(L,x,p)e^-—G(0,x,p) (8) 

For the external regions / and the matching of ip at x = and x = L determine 
the coefficients B and A respectively, so that 

Akx' 



ipAx,p) = iae~ ipx f dx'iaG(x',0,p)e l 

Jo 



" -G(L,0,p)e [kL e- [px ^G(0, 0,p)e~ ip " (9) 



p + k p — k 



and 



$ ni (x,p) =iae ip{x ~ L) [ L dx'iaG(x',L,p)e 

Jo 



ikx' 



G(L, L, p)e ip{x ~ L) e ikL - -^G^, L, p)e ip{x - L) (10) 



p + k ' ' p — k 

m c ip(x-L) c ikL + la 
p2 _ fc2 p2 _ fc2 



Plane wave tunnelling in a square barrier 



5 



For comparison, only the last term in the r.h.s. of (8), the last three ones in (9) and the 
third one in (10) appear in the shutter problem. 

The crucial advantage of the method is that almost all the Green functions involved 
in the equations above can be expanded as a sum of terms which are simple poles in p n . 



3. Analytical structure of the p-dependent solution 

General theorems [I] and the explicit analytical derivation (Appendix B) of G(x,x',p) 
show that \G(x, x',p)\ ~ l/\p\ — > as \p\ — > oo in the complex plane for almost any 
x G [0, L\ and x' G [0, L] . Then the Mittag-Leffler theorem tells that 

G{x,x',p) = Y.^^ 1 (11) 

n P-Pn 

The exception is for G(0, 0,p) and G(L,L,p), the modulus of which grows as \p\ for 
\p\ — > oo in the lower half complex plane p and requires some special care (Appendix 
C). 

The residues in the r.h.s of (11) can be easily computed [3] and one finds 
C n (x,x') = u n (x)u n (x') / N n , where the functions u n (x) belong to the denumerable 
set of the resonant solutions satisfying L x u n (x) = with RBCs 

d x u n (x)\ x=0 = -ip n u n (0) , d x u n (x)\ x=L = \p n U n {L) 

and N n are suitable normalization factors (Appendix A). 
The inverse Laplace transform 

1 /•c+ioc _ \ r+oo . p 2 

^{x,t) = —\ dse st ip(x,p(s)) = / dppe-'amtyfop) , (12) 

Z7T1 Jc-ioo ZTTm J -oo 

written as an integral over the real momentum variable p (Appendix D), leads one to 
consider the pole expansion of p tp(x,p) in the integrand, which can be worked out for 
each region J, // and III (Appendix E): 

P_ u n (0) 

P-Pn N n Jo 

+ a _^ G (L, 0, -k)e~ ipx e ikL - a Y — e " ip V fcL 

P+k n P-Pnk+Pn N n 

+ a?G(0, 0, 0)e~ ipx - J-^—G{Q, 0, k)e~ ipx (13) 
k kp — k 

P 2 1 1 u l(°)„-ipx , ia f 1 , 1 yjh* a -ipxs 



pMx,p) =iae~»"X: P m f dx'u n (x')e ik >' 



+ a y — - - wV y e ~ ipx + — (— — + )(e 

n P-PnPnk-p n N n 2 p - k p+k 



c 



pif) H (x,p) =iaV — - — ^7—/ dx'u n (x')e 
„p-p n N n Jo 



i\ Akx' 



+ a -L-G{L, x, -k)e ikL - a e ikL V Pn - 1 (14) 

p+k V > ^P-Pnk + P n N n K ' 

~ a 7G{0,x,k) +a}_^ 

p-k n P-Pnk-p n N n 



Plane wave tunnelling in a square barrier 



6 



p j m (x, p) = ia eW- z ) £ / L dx' Mn (x') e ifc *' 

nP-Pn N n Jo 

2 

- a V -G(L, L, 0) e ip(x - L) e ikL + -^—^G(L, L, -k)e ip{x ~ L) e ikL 
k p + k k 

_ Q y p2 1 1 u l( L ) c ip(x-D c ikL / 15 n 

n P-PnPnk + p n N n 

n P-Pnk-p n N n k-p 

2 p + k p — k 

In the third row of (14) and in the fourth row of (15) one recognizes the terms of the 
shutter problem given in [2]. The shutter terms in the region I, not calculated in [2], 
are given in the third and fourth rows of (13). 

4. The time-dependent solution 

For each of the terms in the expressions above, the integrals stemming from (12) can be 
brought to the form of an integral representation of the error function erfc(z) so that 
their inverse Laplace transform can be carried out thoroughly (Appendix F). We obtain: 



„ N n Jo 



+ ikG{L,0,-k)e [B_ k \-ie 2,— [B n \ (16) 

n K ' Pn Jv n 

+ iG(0,0,0)± [5j)-iG(o,o,*)i [gfl+iX;^ , * &} 

K K n p n K p n l\ n 

where the factors in square brackets embody the time (r = t / 2m) dependence of the 
solution and are the result of the integrations over the momentum p , namely: 

; ZE 

r ■ EC 2 

A n = -e 1 ^(p n e 2/ -*>» erfc(y^ XjTl ) - 

■a£ 2 

= e^e"- 1 .-* erfc(y- x -k) 

j • 2 

-B n = e^e^— erfc(y- X)n ) 



7TT 



e 1 * x — 2 



56 = — -e* 



7TT 2r 



S£ = - e*(*V*-M erfc(y^ k ) - ^=(k - f )) (17) 

\/7TT zr 



= - e 1 ^ (p£e y -*.» erfc(y„ Xjn ) - ^={p Tl ~ 77-)) 

y 7TT zr 
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S{ = - e~' 1Tk2 erfc{\V^k) 
Si = - e ~ iTk2 erfc(-iV^k) 

s 3 = ~ e^e y -^ erfc(y- Xik ) 

■ x 2 2 

Si = - e^e"— erfc(y_ x „ k ) 

■ 7T 1 

where we have introduced the variables y x>q = e" 1 ^ (4r) _ 5 — 2rq) for 

Q — Pni k, k. 

Likewise 

Mx,t) = -E^#(/ L dx'^)e^')[^] 

+ \kG(L, x, -k)f* m - ie- E u n (L)u n (x) ^ (ig) 

n ft ' Pn -< v n 



where 



where 



i*G(0,s,*) [^]+iE rg; 



A£ J = -p n e- irp: 



2 



rt 



5£i = - e- iTfc2 er/c(-iv / r7A;) 

£ n J = - e~ iTp « erfc{\\/Vrp n ) (19) 

= - e" irfe2 er/c(iv / r7A;) 
S" = -e~ iTp " erfc(iV^p n ) 

Finally 

"71 



„ JO 



iG(L, L, 0)^- [5 m ] + iC(L, L, -fc)^- [5™] 



1C 



1 1 ul{L) 



^EfrrrT^] ( 20 ) 



Pnk+p n N n 

n, (n\n, ( J\ 

nil 



- i*G(0, L, fc) + i £ 

n ft J) n iv n 

- ^ [C[ n ] - U kx [Ci n \ + U kL [Ci n ] + U kL [C{ n \ 



tin 



e 1 *r (p ne yx-L, n erfc(y x - L>n ) + —=) 

\ TXT 



= - e^ffcVU-* erMy.-L.-k) + -^=(-k + ^)) 
K = - e 1 4- (p'e 2 ^.™ erfc(y x - L>n ) + -Wp n + -15—)) 
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gin = _ e i^ e vl. L , k er fc(y x -L,k) 

S 1 " = - e^e^— erfc(y x _ L , n ) (21) 

C[ H = - e" iTfc2 erfc(-iV^k) 

C\ u = - e" iTfc2 erfc(iV^k) 

C\ n = - ^ ij ^ e yl-L,- k erfc(y x - L - k ) 

C{ n = - ^ i£ ^ L e yi-L,H erfc(y x „ L , k ) 

In the equations (16), (18) and (20) the terms [S] are the ones arising in the shutter 
problem. The particular values of the Green function involved are calculated in 
Appendix B. 



5. The short time limit 



The t — > limit is interesting both as a check of the calculation above and for the 
study of the scattered wave at short times. Here we aim only to recover the initial wave 
function ip(x,0) = e lkx at t — 0, which must happen in each of the regions I, II and 
III. 

Careful use of the limiting values of erfc(z) and/or w(z) = e~ z2 erfc(—iz) for z — > 
and for z — > oo in different directions of the complex z plane must be made. Notice 
that y Xtq tends to oo in different directions for different q. Also the properties of the 
set of resonant functions basis of the space of solutions are crucial in region 

II (Appendix G). Asymptotic expressions for erfc(z) and w(z) can be used to obtain 
approximations to the form of ip(x,t) in each region for small values of t, a task that 
will be faced elsewhere. 

In region /, the factor tends to — 2iS(x) and also Sq , S[ and S„ tend to a 
distribution concentrated in x = . The factors B J _ k , B T n , S{ and S{ vanish exactly in 
this limit and both S[ and tend to — 1 . Therefore, as required for x < , 



limt^o ipi(x,t) 



In region II, the factors A^ 1 , B T \ , , S[ T and tend to the value —1 so that 
the [A] term in (18) gives 

[ L dx > V M^K^) = e i kx (22) 
Jo „ N n 

The [B] terms yield 

ikG(L, x, -k)e lkL - ie ikL V P* M^Kfr) 

= ike ikL Y 1 u n (L)u n (x) _ ^ y n u n (L)u n (x) 

^-k-p n N n ^k + p n N n 

= _ ie i fc L^ ^(^K(x) =Q (23) 

n N n 
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and likewise for the [S] terms: 
ikG(0,x,k)-i^—- 

n "> fn 1 ' n 

_ . fe y> 1 U n (0)u n (x) p n U n (0)u n (x) 

= A y U n (0)u n (x) =Q ^ 

n N n 

In region the factor A 1 ^ 1 tends to — 2\8{x — L) and also Bq 11 , B H jf and B n H 
approach distributions concentrated in x = L . Thus they vanish for x > L. The factors 
S IH , S n n , C\ u and C{ H vanish exactly in this limit, whereas both C[ u and C\ u tend 
to —1. Therefore the final result is that, for x > L , also 

limt^o ipin(x,t) = e lhx 

We have quoted singularities at the points x = and x = L respectively in the 
regions / and but these points must be excluded from the domain of respectively 
ipi(x,t) and ipin(x,t) , as argued in Appendix F. 

6. The large time limit 

For r — > oo we see that also y X)Q — > oo as before, but in yet different directions of the 
complex plane for the different q. 

In region I, the factors A l n , B l _ k , B l n , S\ , S l n , S\ and S[ vanish, whereas 

S£ -> - 2k 2 e-' lEkt e~ ikx 

S{ -> - 2e- i£fc * (25) 
S*3 — > — 2e~ LBfct e _ifcr 

Using equation (B.4), some of the e~ lfcr terms cancel out in (16) and one is left with 

lim^oo ^i{x, t) = e - iEkt {e ikx + R(k)e~ ihx ) (26) 
= e~ iEkt <p™(x) 

that is, the asymptotic scattering in solution, where = k 2 / 2m (see Appendix H). 

In region II, the factors A n J , B™ k , B^ 1 and S^ 1 vanish, whereas S^ 1 — > — 2e~ lrfc \ 
Therefore, using (B.2), for < x < L we have 

limt^oo ipn(x, t) = 2ik G(0, x, /c)e" irfc2 

= e-' lEkt (j)l n (x) (27) 

Thus the infinite plane wave evolves into the same final state of the shutter initial 
condition also in region //. 

For the region III, the only non-vanishing factors are 
gin y 2 (j^( x — -k)g — lT ^ 2 

C{ n -> -2e" iTfc2 (28) 
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Then the terms corresponding to C^ 11 and C[ n in (20) cancel each other and, again 
using (B.2), one is left with the same asymptotic solution of the shutter 

limt-^oo i/>iii{x, t) = 2ik G{0, L, k)e ik(x - L) e- wk2 

= T(k)e ikx e~ iEkt (29) 

= e- i£fei C(^) 

The finite time behaviour is thus made up of transient modes which quickly dampen 
out, leaving (for k > 0) the scattering asymptotic solution with outgoing BC at x — L. 

7. Conclusions 

A solution for the time evolution of an infinite plane wave in the presence of a simple 
square barrier has been worked out for each of the regions of the potential and Section 
4 is the main result of this work. Among other terms, this solution contains a sum of 
explicit analytical contributions corresponding to each of the (infinitely many) resonance 
poles. As in previous works in related problems [1] [2], only the location of these poles 
needs to be obtained by numerical methods. 

In the solution obtained the shutter terms have been pinpointed. Similarly, the 
contributions to the time-evolved wave function coming from the segments of the initial 
wave function lying inside and at the right of the barrier can be identified. 

As the main application of this knowledge we envisage the possibility of studying the 
enhancement or the suppression of the transmission of the single Fourier components 
of any realistic wave packet. This should provide new detailed insight on interesting 
phenomena like the super- luminal tunnelling [5], the breakdown of energy conservation 
[6] by transient interference in wave packet collisions with barriers [7j or the rising of 
forerunners. 

The explicit solution obtained for finite time is useful for deriving approximations 
valid for short times, hence for the study of transient structures and forerunners. An 
immediate result of the work is that the large time limit yields the same stationary 
solution of the shutter. This again shows that the resonances contribute only to transient 
structures of the scattered wave. 
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Appendix A. Resonant solutions 

We adopt units such that h = 1 . The homogeneous equation L x u(x) = with the 
RBCs (3) has solutions only for a denumerable set of values p n lying in the lower complex 
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plane p and satisfying the condition 

D{p) = Q 2 e ip ' L - ©V ip ' L = , (A.l) 

where p 1 = \fp 2 — 2mV , © = p + p 1 and = p — p' . One can check that if p n is a 
solution, then — p* n is too, so that these values are in symmetrical locations with respect 
to the imaginary axis. We let the label n take integer values (n ^ 0), with the growing 
positive n indicating the p n with growing real positive part and p- n = —p* n - One finds 
that —j < argp n < and n < argp_„ < ^f. 
The solutions of L x u(x) = are 

^{x) = 6{-x)(-2p' n )e-^ x 

+ 6(x)6(L- x) (©„ e ip '- x - ©„ e~ ip '- x ) (A.2) 

+ 6{x - L) (e n e- ie " L - ©„ e" ie " L ) e iPnX , 

up to an arbitrary multiplicative function of p n . 

The residues C n ( )/N n in (11) correspond to the following choice 

of (complex) "norm" [3] 

N n = i«(0) + u*(L)) + 2p n [ L dx u 2 n (x) (A.3) 

Jo 

which takes the value N n = —8mV(p n L + 2i) for the solutions (A.2). 
Appendix B. Analytical Green function 

The analytical solution to the Green equation L x G(x,y,p) = 5(x — y) for x and y in 
the interval [0, L] , and obeying the RBCs (7), can be directly constructed: 

G(x,y,p) = ^- { 2mV (e^-^+s/)) + e -¥{L-(*+y))) (B.l) 

2p' D(p) 

- ® 2 e- ip ' L (9(y - x)e ip ' iy - x) + 6{x - y)e ip ' {x - y) ) 

- Q 2 e lp ' L {6{y - x)e- ip ' {y - x) + 9(x - y)e- ip ' {x - y) ) } , 



where the symmetry x ^ y is explicit. 

The limit \p\ — ► oo in different directions of the complex plane can be directly read 
out in (B.l). For any values of x and y in the interval [0, L\ the Green function vanishes 
as 1/p or faster in any direction, with the only exception of G(0, 0,p) and G(L, L,p) , 
which grow as p in the lower half plane, while still decreasing as 1/p in the real axis. 

The particular cases G(0, x,p) and G(L, x,p) are related to the scattering solutions 
<p m (x) (Appendix H). From the integral formula J L dx [0 (L X G — 5(x — x')) — G L x <p\ = 
and the use of the BCs for <p and the RBCs for G , one obtains [10] 



G(0, x,p) ^\[^ C(^) (0 < x < L) (B.2) 

G(L, x,p) = — V2tt<P^ 4 n (x) (0<x<L) (B.3) 
2p V m 
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If Imp > , the limits x = in (B.2) and x = L in (B.3) can be taken, namely 

G(0,0,p) = C(0) = y v ^ ( 1 + ^)) = C( L ) = (B.4) 

Appendix C. Pole expansion and substractions 

Whenever the Green function G(x,x',p) vanishes as 1/p in the limit \p\ — > oo, the 
Mittag-Leffler theorem applies and the validity of the pole expansion (11) is assured. 
This is not the case for G(0, 0,p) and G(L, L,p) , and a substraction technique must be 
used This is performed by using the contour integral 

1 f , 1 G(x, x' z) 

— f dz — 

27ri Jr z 2 z — p 

where the contour T = {C Q , C v , C n , C s } consists of small closed paths counterclock-wise 
encircling the poles of the integrand (namely C Q around the double pole at z = 0, C p 
around p and C n around each p n ) plus a large circle C s clock- wise encircling all these 
poles. The factor z~ 2 has been introduced to assure that the integrand goes as z~ 2 even 
for G(0, 0,p) and G(L, L,p) so that the integral over C s vanishes when the radius of the 
circle grows to oo. 

From (C.l) one readily obtains (see also [TT] ) 

G(x, x>, p)= P 2 f: l^^M + [G(a . ^ o) + p d p G(x, P) Ip=o] (C2) 

When the pole expansion (11) holds, (C.2) becomes a trivial identity, whereas for 
G(0, 0,p) and G(L,L,p) a holomorphic (linearly growing with \p\ ) part remains. 



Appendix D. Inverse Laplace transform 

The inverse Laplace transform in the variable s in (12), involves an integration in the 
complex s-plane along a line parallel to the imaginary axis with Re c > 0. In the p-plane 
(recall p = \/iV 2m s) , this path translates to a hyperbola-like one with asymptotes in 

■ v 2 

the positive real and imaginary axis. For t > , the factor e _1 2^* assures that it can be 
deformed into an integration from +oo to — oo along (and slightly above) the real axis 
if the integrand has poles only in the lower half plane (and on the real axis). For t < , 
the path can be closed along a quarter circle of large radius in the 1st quadrant, thus 
enclosing a region without poles and giving ip(x,t) = 0, consistently with the Laplace 
method and causality [4J. 



Appendix E. Pole expansion of pip(x,p) 
We rewrite (8) as 

jnfap) = i« [ L dx / f( I (x' 1 x 1 p)e ikx '-ae ikL fi I (x 1 p)-afi I (x 1 p) (E.l) 
J o 
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where 



vlil , j, CiL.r.p) „ G(0,x,p) 



fi (x',x,p) = G(x', x,p) , f 2 (x,p) = — — and f 3 (x,p) 



p + k p — k 

According to the asymptotic behaviour of the // J (p) (i = 1,2,3) above we consider 
the contour integrals 

= — / dz //'(*) (E.2) 

27ri Jr, z — p 

where the values n x = 0, and n 2 = n 3 = 1 are assigned so that the integrand of (E.2) 
behaves as z~ 2 for large z . The contours Tj include circles around the poles in the 
respective integrand (namely p , p n , k or —k ) plus the large circle C s , as in Appendix 
C. For % — 1, equation (E.2) yields // J (p), which must be multiplied by p later on, 
whereas for % = 2, 3 one directly obtains p// 7 (p). Collecting these results and using the 
pole expansion of the Green function, (14) is obtained. 

The derivation of (13) and (15) follows the same lines. There we define 

ei i \ n( , n x ,// \ G(Q,Q,p) I( , G(L,0,p) 
fi (P) = G(x , 0,p) , f 2 (p) = p _ k , f 3 (p) = , 

and 

//"(„) = G(J,L,p) , Sl-(p) = , fl"(p) = , 

where the powers required re n\ = 0, n 2 = — 1 and 713 = 1 . Notice that the negative 
power n 2 = —1, needed by the asymptotic behaviour of G(0, 0, z) and G(L,L,z), 
introduces an extra pole at z = and is related to the discussion of Appendix C. 

Appendix F. erfc (^-related integrals 

In the computation of ipi(x,t),ipn(x,t) and ipin(x,t), the following types of integrals 



arise: 



: 37T 

f+°° ■ 2 ■ e~ ~ -X 2 

I (x) =— dp e~ 1Tp e lxp = ^=e 1 ^ (F.l) 
m J -00 yvrr 

1; 3tt 

=-/_ oo d PP e-^^=-^-^ (F.2) 

/o(x, q) = - [ + °° dp — e-^e^ = -2 M(x, g, i) = -e^e^- er/c^J (F.3) 
m J-00 p — q 

h(x, q) = - / + °° dp e- Wp2 e ixp = I (x) + ql (x, q) (FA) 

17T J-00 p — q 

/ 2 (x, q) = - dp e~ 1Tp e lxp = h(x) + ql (x) + g 2 / (x, g) , (F.5) 

m J -00 p — q 

— 1 r+°° 1 • t 1 ■ 
M(x, q,t) = — dp e" 1 ^ e lxp 



where 



is the Moshinsky function, often found in the literature. 
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By means of the change of variables p — > p — ^- they can be related to the integral 
representation of the error function [12] 

2 

2 I /-H-oo q— u 

w(z) = e~ z erfc(-iz) = / du Imz > (F.6) 

7T J -co U — Z 

so that 

M(x,q,t) = -e 1 ^ w(ie 4 ]] l^.( x ~ ~ ?)) ■ 

— i 77 

In the expressions above, is convenient the use of the variables y x>q = ^=(x — 2rg) , 
where r = tt~ 

When converting the exponent irp 2 into w 2 , namely u = y/irp , an integration along 
a line running between ±\/ioo results. In the case q = p n , rotating this path by an angle 
7r/4 towards the real axis crosses the pole \f\rp n , contributing a residue. For q = ±k, 
a small negative imaginary part is supposed (see Appendix D). The case Imz < can 
be tackled by using the property w(z*) = w*(—z). 

The results of the integrals above, as given in (17), (19) and (21), are 

A\! = h(0, Pn ) 

A{ n = I x {x-L,p n ) 

B* n = h{x-L) 

Ai = h(-x,p n ) BU£ = I 3 (x-L,-k) 

Bl k = S{ = I (-x,-k) B l n u = I 2 (x-L,p n ) 

Bl n = Sl n = Io(-X,Pn) S{ H = I (x-L,k) . , 

S* = h(-x) Si 11 = I (x-L, Pn ) { • ] 

S{ = I 2 (-x,k) S\ = B u k = C[ n = l {0,-k) 

Si = I (-x,k) S[ = S{ 1 = Ci H = l (0,k) 

C\ n = I (x-L,-k) 

C{ H = I (x-L,k) 

B» = Si 1 = l (0,p n ) 

The computation needs some care in the limit r — > when we simultaneously let x — , 
as is the case for A 1 ^ = 1^0, p n ), because the change of variables p — ► j- becomes ill 

defined. In fact, taking first x — , the naive result A^ = —p n e~ 1TPn erfc(i\nr) — ^7= is 
singular when r — > . The correct result is attainable as follows: The troublesome 
first term in (18) describes the time evolution in region II of the initial condition 
i/jq(x) — e lkx governed by the retarded propagator G R (x,x' ,t) which is the inverse 
Laplace transformed of the G(x,x',p) (subject to RBCs). As such, it has the expansion 

r^R / 1 _iA V"^ —i — t U n (x)u n (x ) 

G K (x,x',t) = '^pn ' , 

hence the result for A™ quoted in (19). Performing the inverse Laplace transform by 
integration in the variable p is not convenient in this case, whereas it is more direct in 

2 

the variable s = = — iE p . 
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The singularities for r — > reported in Section 5 have the same origin and, in this 
sense, the points x = and x = L must be considered as belonging to the non trivial 
region II. 

Appendix G. t — > 0, oo limits of e Vx ><> erfc{y X)q ) 
We see that, for t — ► , 



/m 
2t 



2/a;,g ~e * ^/ — X ^ e * • OO 



regardless of g being equal to ±k, p n or p_ n . 
For t — > oo one has 



where f < n < 2e , -| < 0_ n < , fc = ^ and = -f . 
The following asymptotic formulas [12] apply: 

7T 

Mm z -*oc erfc(z) = < ^ 

lim^oo e z erfc(z) ~ — =— (1 + —^ L ) —> \arq z\ < — 



Only the case 4>k = is out of their range, in which case the relationship 

2 2 2 

e f*,* erfc(y x , k ) = 2e y ^ k - e y ^ k erfc(-y Xtk ) , 
where \arg(-y Xjk )\ = \ , is helpful. 



Appendix H. Scattering solutions 

The usual scattering solutions of L x ip(x,p) = for continuous real energies E p = 
p 2 / 2m > may be characterized according to the homogeneous BC they obey. 
Contrary to the resonant solutions, which are subject to two BCs, each one of the 
scattering solutions solution is subject to one homogeneous BC, alternatively at x = 
or at x = L : 

d x ^\x)\ L =ip<PT(L) 
d x <f>r(x)\o = -ipC(O) 
d x <P° ut (x)\o =lpC(0) 

d x <Pr\*)\L= -ipC(L), (H.1) 



where r(l) labels the right (left)-moving character of the impinging wave. Notice that the 
reverse notation is often found in the literature, with r(l) indicating the wave coming 
from the right (left) of the barrier. 
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These solutions, for instance 

\/27r y p 
are 5-normalized in energy. 
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